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Abstract. Let A C N n be a finite set, and K C R n be a compact semial- 
gebraic set. An A-truncated moment sequence (.4-tms) is a vector y = (y a ) 
indexed by elements in A. The ,4-truncated i^-moment problem (.A-TKMP) 
studies whether or not a given ,4-tms y admits a A"-measure fi, i.e., a nonneg- 
ative Borel measure fi supported in K such that y a = f K x a d/i for all a £ A. 
This paper proposes a numerical algorithm for solving *4-TKMPs. It is based 
on finding a flat extension of y by solving a hierarchy of semidefinite relax- 
ations {(SDRjfcj^Lj^ for a moment optimization problem, whose objective R is 
generated in a certain randomized way. If y admits no i^-measures and IR[:r]_4 
is K-iull (there exists p = "^2 aej ^PaX a that is positive on K), then (SDR)j, is 
infeasible for all k big enough, which gives a certificate for the nonexistence 
of representing measures. If y admits a i\-measure, then for almost all gener- 
ated R, we prove that: i) we can asymptotically get a flat extension of y by 
solving the hierarchy {(SDR^jJ^Lj ; ii) under a general condition that is almost 
sufficient and necessary, we can get a fiat extension of y by solving (SDR)fc 
for some fe; this occurred in all our numerical experiments; iii) the obtained 
flat extensions admit a r-atomic i^-measure with r < \A\. The decomposition 
problems for completely positive matrices and sums of even powers of real lin- 
ear forms, and the standard truncated i^-moment problems, are special cases 
of .A-TKMPs. They can be solved numerically by this algorithm. 



Let A C N" be a finite set (N is the set of nonnegative integers). An A-truncated 
moment sequence (A-tms) is a vector y — (y a )a£A in K."^ (the space of real vectors 
indexed by elements in A). For a :— {a\, . . . ,a n ) € N ra , denote \a\ := oti + ■ ■ ■ + a n . 
The degree of A is deg(.4) := max{|a| : a E A}. Let K be the semialgebraic set 



defined by two tuples of polynomials h — (h\, . . . , h mi ) and g = (gi, . . . ,g m2 )- A 
nonnegative Borel measure fi on R™ is called a K -measure if its support, denoted by 
supp(/i), is contained in K. For x := (xi, . . . , x n ) G E n and a := (ai, . . . , a n ) G N n , 
denote x a :— x" 1 ■ ■ ■ x"'* . An ,4-tms y is said to admit a X-measure /x if 



The measure /x satisfying the above is called a K -representing measure for y. Let 
measly, K) denote the set of all X-measures admitted by y. Denote 
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1. Introduction 






M A (K) := {y <E R : meas(y, K) ^ 0}. 
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The ^4-truncated if -moment problem (^4-TKMP) is to determine whether a given 
A-tms y admits a if-measure or not? If it does not, can we get a certificate for 
that? If it does, how can we obtain a if-representing measure? Preferably, we are 
often interested in finitely atomic ^-representing measures. (A measure is finitely 
atomic if its support is a finite set, and is r-atomic if its support consists of at 
most r distinct points.) This paper is to present a numerical algorithm for solving 
ATKMPs. 

1.1. Two special cases. Many interesting computational problems can be formu- 
lated as .A-TKMPs with appropriate A and K. Here we list two of them. 

The first one is the decomposition problem for completely positive matrices 
(cf. [5]). A symmetric n x n matrix C is called completely positive if there ex- 
ist vectors u\, . . . ,u r € R" (the nonnegative orthant of K n ) such that 



The above is called a CP-decomposition of C, if it exists. How can we determine 
whether a matrix C is completely positive or not? If it is not, can we get a certificate 
for that? If it is, how can we get a CP-decomposition for CI As we will see in 
Section [pj this problem can be formulated as an ,4-TKMP with 



The second one is the decomposition problem for sums of even powers (SOEP) 
of real linear forms (cf. 35 ). (A form is a homogeneous polynomial.) A form / of 
an even degree m is SOEP if there exist real linear forms L%, . . . , L r such that 



The above is called an SOEP- decomposition of /, if it exists. How can we determine 
whether a form / is SOEP or not? If it is not, can we get a certificate for that? If 
it is, how can we get an SOEP-decomposition for /? As we will see in Section [51 
this problem can also be formulated as an „4-TKMP with 

A = {a e N n : \a\ = to}, K = {x e E n : x T x = 1, xi + • • • + x n > 0}. 

It is typically quite difficult to detect the existence of or computing CP/SOEP- 
decompositions. In the prior existing work, there are no much efficient numerical 
methods for solving such decomposition problems (except some special cases), in 
the author's best knowledge. In this paper, we will show how these problems can 
be solved numerically as special cases of „4-TKMPs. 

1.2. Standard truncated if-moment problems. Denote NJJ := {a € N" : 

M < d}. When A — NJJ, the ,4-TKMP is specialized to the standard truncated 
JC-moment problem (TKMP), and an N^-tms is called a tms of degree d. Curto 
and Fialkow originally studied TKMPs and have made foundational work in the 
field. We refer to [U [TTJ [12] and the references therein. Here we give a short review 
for TKMP. 

Every tms z € M Nd defines a Riesz functional ££ z acting on M.[x]d (the space of 
real polynomials in x := (x±, . . . , x n ) of degrees at most d) as 



C = uiuj H + u r uj . 



A = {a e N n : \a\ = 2}, K = {x e R+ : x x + ■ ■ ■ + x 



!}• 



(1.2) 
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We also denote (p,z) := ^f z (p). The Riesz functional Ji? z is A -positive if 

J?z{p)>0 Vp6%] d : p\ K >0. 

The A-positivity of Jz? z is necessary for z to admit a A-measure. When K is 
compact, it is also sufficient, which can be implied from the proof of Tchakaloff 's 
Theorem [39]. However, it is typically very difficult to check whether a Riesz 
functional is A-positive or not. 

A more favorable condition than A"-positivity of J£ z is flatness. For convenience 
of description, suppose A — W 1 . We denote X > (resp., X >- 0) if the matrix 
X is symmetric positive semidefinite (resp., definite). For a tms z £ R N 2fc, define 
Mk{y) to be the symmetric matrix, which is linear in z, such that 

^z{p 2 ) = p T M k (z)p Vp £ R[x] k . 

(For convenience, we also use p to denote the vector of coefficients of p(x) in the 
graded lexicographical ordering.) The matrix M k (z) is called a k-th order moment 
matrix. If z admits a A-measure /i, then 

p T M fe (z)p = ^ 2 (p 2 ) = / p 2 dfi>0 V P ER[x} k . 

J K 

This implies that 

(1.3) M fe (z) b 0. 

Hence, (|1.3[) is necessary for z to admit a measure on R™, but typically not sufficient. 
However, if (|1.3p is satisfied and z is flat, i.e., 

(1.4) rankMjt_i(z) =rankM fe (^), 

then z admits a unique measure, which is r-atomic with r = rank Mk (z) . When 
A is a semialgebraic set as in (ll.l[) . there is a similar version of this result (cf. 
Theorem I2.2|) . This is an important result of Curto and Fialkow (cf. 11J). For 
convenience of notion, when A = K™, we simply say z is flat if it satisfies ()1.4j) and 
IjOj) . (When A is as in (L~I]> . we say z is flat if it satisfies {13]) and (pT2]l .) 

Flatness is very useful for solving truncated moment problems. Let y £ M.^ d , z € 
IR N '! be two tms'. We say y is a truncation of z, or equivalently, z is an extension of 
y, if d < e and y Q = z a for all a £ NJJ. We denote by z\d the subvector of z, whose 
entries are indexed by a € NJJ. So, z is an extension of y if and only if y = z\d- 
If z is flat and extends y, we say z is a flat extension of ?/. Similarly, if w is an 
extension of z and z is flat, then we say z is a flat truncation of w. Clearly, if y is 
not flat but has a flat extension w, then y and w commonly admit a finitely atomic 
measure (cf. Theorem 12. 2[) . Indeed, Curto and Fialkow [TT] further proved that: a 
tms y £ E N S admits a A-measure if and only if it has a flat extension w £ K N 2fc 
(cf. Theorem 12. 3p . However, little is known on how to get flat extensions. 

When A is compact as in (jl.ip . Helton and Nie [20] proposed a semidefinite 
approach for solving TKMPs. Its basic idea is to find flat extensions through 
semidefinite relaxations with moment and localizing matrices. They proved that: 
if y £ admits no A-measures, then a semidefinite program will be infeasible; if 
y £ M^ d admits a A-measure and there exists fj, £ meas(y, A) such that supp(/x) C 
Z(p) (Z(p) denotes the real zero set of p), where p is a polynomial having certain 
weighted SOS type representations and Z(p) is finite, then a flat extension of y can 
be found by solving a semidefinite program. However, in more general cases, no 
much were known on how to get flat extensions. 
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1.3. Contributions. This paper studies the new and broader class of moment 
problems: ^4-TKMPs. The CP/SOEP-decomposition problems and standard trun- 
cated if-moment problems are special cases of „4-TKMPs. 

For an .A-tms y, a tms w £ M. Nzh is an extension of y (or y is a truncation of w) 
if w a — y a for all a £ A. Denote by w\a the subvector of w whose indices belong 
to A. Clearly, if w\a = y and w is flat, then w and y commonly admit a finitely 
atomic measure. When does y admit a i^-measure? If it does, how can we get such 
a measure? If it does not, can we get a certificate for the nonexistence of such a 
measure? These are the main questions in ,4-TKMPs. 

This paper proposes a numerical algorithm (i.e., Algorithm I4.2[) for solving A- 
TKMPs. It is based on solving a hierarchy of semidefinite relaxations (we denote 
by {(SDR)^}^^ here for convenience), for a moment optimization problem whose 
objective is a Riesz functional (R,w) :— J*? W (R). The objective R is generated in 
a certain randomized way. Assume K is given as in (jl.ll) and is compact. Denote 
R[£]*4 : = span{a; Q : a £ A}. If y admits no X-measures and R[:e].a is Jf-full (i.e., 
there exists p £ M.[x]a that is positive on K), then (SDR)fc will be infeasible for 
all k big enough. This gives a certificate for the nonexistence of a i-T-representing 
measure for y. If y admits a i^-measure, then, for almost all generated R, this 
algorithm has the following properties: 

i) We can asymptotically get a flat extension of y by solving the hierarchy 
{(SDR)fc}^ 1 . So, the convergence is guaranteed with probability one. 

ii) We can get a flat extension of y by solving (SDR)^ for some k, under a 
general condition that is almost sufficient and necessary. This implies that 
the finite convergence is very likely to happen. 

iii) The obtained flat extensions admit a r-atomic if-measure with r < \A\. 

In all our numerical experiments, the finite convergence was always observed, and 
we got r-atomic if-representing measures with r < \A\ when they exist. 

CP/SOEP-decomposition problems are special cases of ,4-TKMPs. So, this al- 
gorithm can be applied to solve them. If such decompositions do not exist, then 
the resulting (SDR)^ will be infeasible for all big k, which gives a certificate for the 
nonexistence; if they exist, we can get such decompositions, either asymptotically 
or in finitely many steps (very likely to happen). This algorithm can also solve stan- 
dard TKMPs. In the author's best knowledge, this is the first numerical algorithm 
that has the aforementioned properties. 

The paper is organized as follows. Section [5] reviews some backgrounds; Section[3] 
presents some properties of ,4-TKMPs; Section @] describes the algorithm; Section[5] 
proves the convergence results; Section [6] gives applications in solving CP/SOEP- 
decomposition problems and the standard TKMPs. 

2. Backgrounds 

Notation The symbol N (resp., R) denotes the set of nonnegative integers (resp., 
real numbers). For ( £ K, \t] (resp., [t\) denotes the smallest integer not smaller 
(resp., the largest integer not greater) than t. For denote by [x]a = {x a ) a ^A 

the vector of monomials, whose exponents are from A, ordered in the lexicographical 
ordering. Denote N% ■= {a £ W 1 : \a\ < d}. When A = the vector [x]a is 
denoted by [x]d- Denote [k] := {1, . . . , k}. The symbol R[x] — R[xi, . . . , x n ] denotes 
the ring of polynomials in x := (x\, . . . , x n ) with real coefficients. 
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For a set S C R n , \S\ denotes its cardinality, and int(S) denotes its interior. 
The symbol Vd(S) denotes the set of polynomials in WL[x]d that are nonnegative 
on S. The superscript T denotes the transpose of a matrix. For u £ M. N , denote 
||w||2 := V u T u and B{u 1 r) :— {x € R™ : ||x — u|| 2 < r}. For a polynomial p <E R[x], 
||p||2 denotes the 2-norm of the coefficient vector of p. Denote by S fc = {x £ R fe+1 : 
ll^lb = 1} the fc-dimensional unit sphere. Denote by S N the space of N x iV real 
symmetric matrices. For a matrix A, \\A\\2 denotes its standard operator 2-norm, 
and ||^4||f denotes the Frobenius norm of A, i.e., \\A\\p — y/Trace(A T A). 

2.1. Standard truncated if-moment problems. Let z £ R^ 1 and K C R™. 

Bayer and Teichmann [T] proved that: z admits a if-measure if and only if it 
admits a r-atomic if -measure v with r < ■ A nice exposition for this result can 

be found in Laurent [27J Theorem 5.8]. When K is compact, we can characterize 
the existence of representing measures via Riesz functionals. This can be implied 
by the proof of Tchakaloff's Theorem [39] , 

Theorem 2.1 (Tchakaloff ) . Let K be a compact set in W 1 . A tms z £ M Nd admits 
a K -measure if and only if its Riesz functional J/f 2 is K -positive. 

When K is noncompact, the above might not be true. A stronger condition is 
that Jzfj, is strictly K -positive, i.e., 

J? y (p)>0 Vp£R[x} d :p\ K >0,p\ K ^0. 

When K is a determining set (i.e., p = whenever p\x = 0), if ££ z is strictly K- 
positive, then z admits a if- measure (cf. [171 Theorem 1.3]). Typically, checking 
if-positivity or strict if-positivity is quite difficult. 

A more useful condition than if-positivity is the positive semidefiniteness of 
localizing matrices. For z £ R N2fc and h £ R[x]2fe, define L q (z) to be the symmetric 
matrix, which is linear in z, such that 

(2.1) J? z (qp 2 ) = p T (Lf\z))p Vp€R[ar] fc _ rdog(g)/21 . 

The matrix L q k \z) is called the fc-th localizing matrix of q generated by z. When 
5=1, Lq k \z) coincides with the moment matrix Mk(z). 

Let K be as in (jl.ip and go = 1. If z £ R Nafc admits a if -measure then 

(2.2) L^(z) = 0(i = l,...,mi), L^(z)h0(j = 0,l,..-,m 2 ). 
This is because 

p T L {k) (z)p = J? z (h i p 2 ) = [ h lP 2 d f i = 0, 

JK 

p T L<£>(z)p = J? z {g oP 2 ) = g 3 p 2 ^ > 0, 

for all p £ M.[x] with m&Xij{deg(hip 2 ),deg(gjp 2 )} < 2k. Therefore, (|2.2|) is neces- 
sary for z to admit a if -measure. Typically, it is not sufficient. Let 

(2.3) d K := max {1, \deg(hi)/2], \deg( 9j )/2]}. 

If, in addition to (|2.2[) . z satisfies the rank condition 

(2.4) rankA4_ dK (z) = rankAf fe (z), 



6 



JIAWANG NIE 



then z admits a unique if-measure, which is finitely atomic. This is an important 
result of Curto and Fialkow. For convenience of notion, we simply say z is flat if z 
satisfies (|2T2"j) and (|2l| . 

Theorem 2.2 ([11]). Let K be defined in ifO]) . If atmsze R N ?<= is flat (i.e., p^] ) 
and i^. £)) hold), then z admits a unique K-measure, which is rank Mk(z)- atomic. 

Flatness is very useful for solving truncated moment problems, as shown by 
Curto and Fialkow [9] [10] . A nice exposition for this can also be found in Laurent 
[26] , For a flat tms, its finitely atomic representing measure can be found by solving 
some eigenvalue problems, as shown by Henrion and Lasserre |21j . 

Clearly, if a tms y is not flat but has a flat extension z, then y admits a K- 
measure. Thus, the existence of a if-representing measure for y can be determined 
by investigating whether y has a flat extension or not. Indeed, Curto and Fialkow 
[IT] proved the following important result. 

Theorem 2.3 ([II]). Let K be as in LP]) . Then a tms y G R N S admits a K- 
measure if and only if it has a flat extension z 6 M N2fc . 

There are other necessary conditions for admitting if -representing measures, like 
the recursively generated relation. We refer to Curto and Fialkow 8, 91 110). 

2.2. Ideals, varieties and positive polynomials. A subset / C M.[x] is called an 
ideal if I + I C / and p-q € L for all p € I and g £ R[4 For a tuple p = (j>±, . . . ,p m ) 
of polynomials in M[x], denote by I(j>) the ideal generated by pi, . . . ,p m , which is 
the set piM[a;] + • • ■ +p m M[a;]. A variety is a subset of C n that consists of common 
zeros of a set of polynomials. A real variety is the intersection of a variety with M. n . 
Given a tuple p = (p\, . . . ,p m ), denote 



A polynomial / G R[x] is a sum of squares (SOS) if there exist fi,...,fk € M.[x] 
such that f = fi + • ■ ■ + fk- The set of all SOS polynomials in n variables and of 
degree d is denoted by £ n .<j. It is a convex cone in IR Nd and has nonempty interior 
for any even d > 0. We refer to Reznick [34] for a survey on SOS polynomials. 

Let K be as in (jl.ip . defined by the polynomial tuples h = (h\, ■ ■ ■ , h mi ) and 
g = (.gi, ■ ■ ■ , g m2 )- Denote g = 1 and 



Clearly, 1(h) is the union of all l2k(h), and -?2fc(M is a truncation of 1(h) with 
degree 2k. The union of all sets Qk{g), denoted by Q(g), is called the quadratic 
module generated by g, and each Qk(g) is a truncation of Q(g) with degree 2k. 

Clearly, if / € 1(h) + Q(g), then / is nonnegative on K. The archimedean 
condition for (h,g) is that there exists N > such that N — ||x||| G 1(h) + Q(g). 
When K is compact, the archimedean condition can be forced to be held by adding 
a redundant condition like N — \\x\W > to the description of K . 

Theorem 2.4 (Putinar, [29p. Let K be as in Suppose the archimedean 

condition holds for (h,g). If f € M[x] is positive on K, then f € 1(h) + Q(g). 



Vc(p) := {v e C" : p(v) = 0}, V R (p) := {v e M n : p(v) = 0}. 



(2.5) 




(2.6) 
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2.3. Semidefinite optimization with moment variables. Scmidcfinite pro- 
gramming (SDP) is very useful in solving moment problems. We refer to [40] 
for a survey on SDP, and refer to [251 1271 120] for semidefinite programs arising from 
moment problems. 

Let be the cone of positive semidefinite matrices in S N . Let JC = x • • • x 
S+ ( be a cone in the space S Nl x • • • x S Ne . A tuple X — (Xi, . . . , Xg) belongs to 
JC if and only if each X, G . A general semidefinite program is 

(2.7) min Trace(CX) s.t. F{X) = f, X G K 

with C a constant tuple in S Nl x • • • x S Nl , J- a linear operator from S Nl x • • • x S Ne 
to a vector space K m , and / a given vector in K m . Let T* be the adjoint operator 
of T . The dual optimization problem of (|2.7[) is 

(2.8) max fy s.t. F*(y) + Z = C, Z G JC. 

Let h and g be the tuples of polynomials describing K as in (fT!]) . Denote 

(2.9) ^ k (g):={w€R^\L^(w)hO, J = 0, 1, . . . , m 2 } , 

(2.10) E k (h) := {w G M N " L$(w) = 0, i = l,...,mi}. 

It can be shown that $&(<?) is the dual cone of Qk(g) and Ek(h) is the dual cone of 
/afc(ft) (cf. [251 [27]). Indeed, E k {h) is a subspace of R N ^. 
A typical semidefinite program in this paper is 

(2.11) min c T w s.t. w\a = y, w G $k(g) n Ek(h) 

w 

with given c G K. N ' i * (deg(^4) < d < 2k) and y G K" 4 - The dual optimization problem 
of (12~TTT) is 

(2.12) max {p,y} s.t. c- p e I 2k (h) + Q k (g). 
P m[x] A 

Any objective value of a feasible solution of (|2.11[) (resp., (|2.12l) ) is an upper bound 
(resp., lower bound) for the optimal value of the other one (this is called weak 
duality). If one of them has an interior point (for (|2.12| it means that there exists 
p G R[a;].4 such that c — p lies in the interior of (Qk{g) + l2k(h)) H and for 

(|2.11| it means that there is a feasible w with Lg){w) >~ for all j), then the other 
one has an optimizer and they have the same optimal value (this is called strong 
duality). We refer to [2] Section 2.4] for duality theory in linear conic optimization. 

3. Properties of ,4-TKMPs 

This section presents some properties of ,4-truncated K-momeiit problems. They 
will be used in the latter sections. 

3.1. ,4-Riesz functionals. An ,4-tms y defines an ,4-Riesz functional Jzfj, that acts 
on R[x].4 as 

aeA aeA 

The vA-Riesz functional J£ y is said to be K -positive if 

J? v (p)>0 Vp G R[x] A : p\ K > 0, 
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and strictly if-positive if 



J2?„(p) > Vp e R[x] A : p\ K > 0, p\ K # 0. 



Clearly, if y admits a if-measure /x, then Jz? y must be if -positive because 



whenever p 6 R[x].4 and £>|k > 0. So, the if-positivity of Jzfj, is a necessary 
condition for y to admit a JT-measure. Indeed, it is also a sufficient if K is compact 
and R[x].a is K-full (i.e., there exists p E M.[x]a such that p > on if). This is a 
result of Fialkow and Nie [18] . 

Theorem 3.1. ( ^TE\ Theorem 2.2]) Suppose K C R™ is compact and M\x\a is 
K-full. Let y be an A-tms such that Jz? y is K -positive. Then, there exist N < 
dimR[x].4 ; Mi, . . . , un G K , and ci, . . . , Cn > 0, such that 



The above theorem immediately implies the following. 

Corollary 3.2. Suppose K C l n is compact and M[x\ A * s K-full. If an A-tms y 
admits no K-measures, then there exists p g R[x].4 such that 



There is a similar version of Theorem 12.31 for ,4-TKMPs. Recall that a tms 
w G R N ™ is called fiat if it satisfies (|2J2]) and (j2T4|t . 

Proposition 3.3. (i) Let K C R" &e a set. Then, an A-tms y admits a K -measure 
if and only if it admits a r -atomic K -measure with r < \A\. 

(ii) Let K be as in Then, an A-tms y admits a K-measure if and only if y 

is extendable to a flat tms w € R^ 2 * for some k. 

Proof, (i) The "if direction is obvious. The "only if direction can be proved by 
a formal repetition of the proof of Theorem 5.8 in Laurent [2 7) . The inequality 
r < |„4| is implied by Caratheodory's Theorem. Here we omit it for cleanness of 
presentation. 

(ii) "=>" Suppose y admits a measure on K. By item (i), it admits a r-atomic 
if-measure, say, y — ci[ui]a + • • • + c r [u r ]j[ with all Ui G K, > and r < \A\. 
For k big enough, the tms w = ci[wi]2fe + ■ • • + c r [u r ]2k € R N 2fc is fiat. Clearly, w is 
an extension of y. 

"<=" Suppose w € R N2fc is flat and w\a = y- By Theorem 12.21 w admits a 
if-measure [i and w a — J K x a d/i for all a € N^- Since w\a = y, we have y a = 
J K x a dfi for all a £ A, i.e., y admits the if- measure fx. □ 

3.2. Extremal extensions. For d > deg(^4), define the d-th extension of y as 



Clearly, £d{y, K) is convex. If meas(y, K) ^ 0, then £d{y, K) ^ 0. 

Lemma 3.4. Let y e R" 4 . If K is compact and R[x]^ is K-full, then £d(y,K) is 
a compact convex set. 




V = ci[uiU H h cat[uat]^. 



Pk > 0, J2^(p) < 0. 



(3-1) 



£ d (y,K) = {z£& N n(K): z\ A = y}. 
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Proof. The compactness of K implies (K) is closed, which can be implied by 
Theorem l2.ll So, £d(y, K) is also closed. If meas(y, K) = 0, then £d(y, K) = and 
we are done. Thus, we can assume meas(y, K) ^ and £d(y,K) ^ 0. We need 
to prove that £d(y,K) is bounded. Since M[x]^4 is JT-full and K is compact, there 
exist p £ R[x].a and e such that p\x > e > 0. There exists M > such that for all 
x £ K 

-M <x a <M Va G NJ. 
For all z € £d(t/, if), the if-positivity of J£ z implies that 

-Mz <z a < Mz Va G N£. 
Similarly, > e implies ^f y (p) > ezo- So, 

|*d < MJSf„(p)/e VaeN^. 
Hence, £d(y,K) is bounded, which completes the proof. □ 

Lemma 3.5. Lei if C R n , d > deg(_4) and y G &j\[[K)- Suppose z is an extreme 
point of £d(y,K). If fi £ meas{z,K), then /j. must be r-atomic with r < \A\. 

Proof. Choose an arbitrary fi £ meas(z, K). If /i is r-atomic and r < \A\ 1 then we 
are done. We derive a contradiction if either /i is r-atomic with r > \A\, or /x is not 
finitely atomic. 

First, consider the case that /i is r-atomic with r > \A\, say, 

/i = ci5(ui) + ■ ■ • + c r 5(u r ), 

with distinct points ui, . . . , u r G K and ci, . . . , c r > 0. Here S(v) denote the Dirac 
measure supported on the point v. For each a £ NJJ, denote 

tuW := «,..., <) T e R' r . 
We show that there must exist f3 £ N2\A such that the system 

(3.2) (w^Y t = 0(Va G A), (w^ t^O 

has a solution t — (ti, . . . , t r ). Suppose otherwise, then for all [3 £ N^\A, 

(w (a) Y t ^ (y a £ A) => (w^Yt = 0. 

This implies that, for every /3 G the vector is a linear combination of 

the vectors w^ a \a £ A), i.e., there exist real numbers pp, a such that 

- (/3) = £^ (Q) - 

For each (3 £ N^\A, let 

Pf3 ■= X - ^ P/3,aX a . 

Note that pp(ui) = for all i — 1, . . . , r. Let J be the ideal generated by pp(0 £ 
N^\A). Then J is zero-dimensional, and the dimension 

D := dimM[a;]/J < dim K [a;] .a = 

By Proposition 2.1 of Sturmfels [38], the number of common zeros of polynomials 
P/3 (/3 £ Nj\A) is equal to D, counting multiplicities. Hence, the cardinality of 
Vc(J) is at most |.4|. However, the distinct points u\,...,u r all belong to Vc(J) 
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and r > D, which is a contradiction. So, there exists t — (ti, . . . ,t r ) satisfying 
(|3.2p . Then, (ci, . . . , c r ) ± e(ii, . . . , t r ) > for e > sufficiently small. Let 



i=l i=l 

They are all nonnegative Borel measures supported in K. Let 



zi = / [^JddMij z 2 = / Fj d d^ 2 - 



Then, both zi and z 2 belong to £d(y,K) and z = \{z\ + 22)- The inequality in 
implies that z\ ^ z, z 2 7^ 2. This contradicts that z is extreme in £d(y, K). 



Second, consider the case that /i is not finitely atomic. Then |supp(/i)| = +00. 
Choose \A\ + 1 distinct points, say, v\, . . . , U1.41+1, from supp(/x). The support of /1 
is the smallest closed set S such that fj,(R n \S) = (cf. [10 Section 4]). So, there 
exists e > such that fi(B(vi, e)) > for all i and the balls B(y\, e), . . . , B{v\a\ + \, e) 
are disjoint from each other. Let Ti :— B(vi, e) n supp(/j,) for i = 1, . . . , \A\, and 

1-41 

T \A\+i '■= supp(^) - (J B(vi,e). 

1=1 

This results in the decomposition supp(^) = T\ U • • • U T\^\ +1 . Note that £t(7i) > 
0, . . . , /i(T|^| +1 ) > and T; fl 2} = whenever i 7^ j. For each j = 1, . . . , \A\ + 1, 
let Hj = (J-Itj, the restriction of /i on Tj. Then 



/ [x] d dm H / [x] d dfi\ A \ +1 . 

Jk Jk 



Each tms J K [x]dd/ii admits a finitely atomic measure supported in Ti (cf. |27l 
Corollary 5.9]). Hence, there exists a measure 9 € meas(z, K) that is r-atomic 
with r > \A\. Therefore, a contradiction can be obtained as in the first case. 

Combining the above two cases, we know the conclusion is true. □ 

3.3. Linear optimization over £d(y,K). Let d > deg(.4) and R £ M.[x]d- Con- 
sider the linear moment optimization problem 

(3.3) min (R,z) s.t. z\ A = y, z e StAK). 

z 

The feasible set of fl3.3[) is £d(y,K). Clearly, if z* is the unique minimizer of (|3.3p . 
then z* is an extreme point of £d{y 7 K). Recall that Vd(K) = {p S M.[x]d ■ p\k > 0}. 
The dual optimization problem of p. 31) is 

(3.4) max (p,y) s.t. R- P eV d (K). 



This is because the dual cone of Md{K) is Pd(K), for any compact set K. This was 
shown by Tchakaloff [39] (also see Laurent [UJ Section 5.2]). 

Proposition 3.6. Let K C R™ be a compact set and y be an A-tms in & A (K). 

(i) // either R\k > or M.[x] A is K-full, then US. 3\) and {3.J$ have the same 
optimal value and iS. 3]) has a minimizer. 

(ii) If J? y is strictly K -positive, then \3.$ has a maximizer. 
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Proof, (i) If R\k > 0, then the origin is an interior point of (|3.4j) . If M[a;]^ is K-fuH, 
then there exists q G R[:e].a C M.[x]d that is strictly positive on K, so (|3.4[) also has 
an interior point. The problem (|3.3[) is feasible, because y G ^^(if). Thus, the 
strong duality holds and the conclusion is true (cf. [SJ Section 2.4]). 

(ii) From y G Mj^K), we know y has a flat extension u> 6 R n 2j» for some A; 
(cf. Proposition [33]). Then the truncation w\d is a feasible point for (|3.3I) . By weak 
duality, the optimal value of (|3.4[) is finite, say, 77. Clearly, the feasible set of Q3.4p 
is closed. Let {pk} C K[x]^ be a sequence such that each R— pk € Vd{K) and 

(Pfc; y) ~* Vi as fc — > 00. 

Let Si = {/ € R[a;].A : J\k = 0}, and S2 be the orthogonal complement of Si in 
R[x] A - So, R[x] a = Si + S 2 , Si _L S 2 . Because (/,y) = and /| K ee for all 
/ G Si, the pfc's in the above can be chosen as pk G S2 and R — pk G Vd(K). 

If the sequence {p&} is bounded, then any of its accumulation points is a max- 
imizer of (|3 .4[) and we are done. Suppose otherwise {pk} is unbounded, say, 
\\PkW2 00 < Let p/j = p/c/Hpfclb be the normalization. The sequence {fib} is 
bounded, and we can generally assume pk — > p* G S 2 . Clearly, ||p*||2 = 1 and 

(Pk,y) = (pk,y)/\\pkh ->■ 0, as £; — s- 00. 

This implies that (p*,y) = 0. From (iJ - p fc )/||p*||2 G V d (K), we get -p*\ K > 0. 
Since p* G S 2 and ||p*||2 = 1, we know ~P*\k ^ 0. The strict i(T-positivity of jSf w 
implies {—p*,y) > 0, which is a contradiction. Thus, the sequence {pfe} must be 
bounded, and the proof is complete. □ 

4. A SEMIDEFINITE ALGORITHM FOR ,4-TKMPS 

In this section, we present a numerical algorithm for solving ,4-truncated K- 
moment problems. To determine whether an .A-tms y admits a ii'-measure or not, 
by Proposition 13.31 it is equivalent to investigating whether y has a flat extension 
or not. If it does not exist, then y does not admit a if-measure. If it exists, then 
we can get a finitely atomic representing measure for y through the flat extension. 

The extension set £d(y,K), defined in (|3.1I) . is very useful in getting flat exten- 
sions. By Lemma 13.51 extreme points of £d{y, K) admit only r-atomic if-measures 
with r < |„4|. Clearly, if (|3.3p has a unique minimizer z* , then z* is an extreme 
point of £d{y, K). A very useful fact is that when we minimize a generic linear func- 
tion over a compact convex set, the minimizer is unique (cf. |36[ Theorem 2.2.4]). 
If K is compact and M[a;]yv is JT-full, the set £d(y, K) is compact convex by Propo- 
sition [33J and (|3.3j) has a minimizer for all R. If K[a;]^ is not if -full, we typically 
need to choose R that is positive definite on K, to guarantee (I3.3[) has a minimizer. 
Therefore, to get an extreme point of £d(y,K), it is enough to solve (|3.3[) for a 
generic positive definite R, no matter R[x]^ is A'-full or not. 

The cone M&{K) is typically quite difficult to describe, and generally we can not 
solve p.3p directly. Recently, there is much work on solving moment optimization 
problems like J33]) by semidefmite relaxations. We refer to [H [201 HI 113 • The 
basic idea is to approximate the cone Md(K) by semidefmite programs. So, we use 
semidefmite relaxations to solve (|3.3I) . This produces a semidefmite algorithm for 
solving ATKMPs. 

Suppose K is compact and K C B(0,p). Denote b := p 2 — WxW 2 ,. Recall that 
h = (hi, . . . , h mi ) and g — (gi, . . . , g m2 ) are the tuples of polynomials describing K 
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as in fjl . 1 [) . For convenience, denote 

9b ■= (91, ■ ■ -,g m2 ,b). 

The set K can be equivalently described as h(x) = 0,gB (x) > 0. Recall the 
definitions of Qk(g) in (|2.6j) and its dual cone $k(g) in (|2.9j) . Note that 

Qk(g B ) = Qk(g) + Q k (b), <s>k(g B ) = n $*(&). 

Let k > d/2 be an integer. The fc-th order scmidchnitc relaxation of f|3 . 3[) is 

(4.1) min (R,w) s.t. w\ A = y, w £ $ fe (ss) n E k (h). 

w 

The dual optimization problem of (|4.ip is 

(4.2) max (p,y) s.t. R- p G I 2k {h) + Qkigs)- 
P eR[x] A 

For every w feasible for (|4.ip and every p feasible for (|4.2[) , we have 

(R,w) > (p,y), 

by weak duality. Thus, the optimal value of (|4.1[) is always greater than or equal 
to that of (|4T2|) . 

Proposition 4.1. Let y £ R A , d > deg(^4), and K C B(0, p) be as inU.l\). 

(i) If w* is a minimizer of and has a flat truncation w*\2t with It > 
deg(^4), then y admits a finitely atomic K -measure \i. 

(ii) If is infeasible for some k, then y admits no K-measures. 

Proof, (i) Let z = w*\2t- Then z is a flat tms. So, z admits a finitely atomic 
JC-measure fi, by Theorem 12.21 Since 2t > deg(^4), z\a — y and z is an extension 
of y. Thus, y also admits the measure fi. 

(ii) Suppose otherwise y admits a ii'-measure. By Proposition 13.31 y can be 
extended to a flat tms w\ £ R N2t: i satisfying (|2.2p . By Theorem 12.21 the tms w\ 
admits a r-atomic X-measure, say, wi — ci[ui]2fc 1 + ■ ■ ■ + c r [u r ]2fc! with all c% > 
and Ui £ K. Let w = c\[u{\2k + ■ ■ ■ + c r [u r ]2k, then w is feasible for (|4.1| . a 
contradiction. So, y admits no X-measures. □ 

Proposition 14.11 can be applied to determine whether an .A-tms y admits a K- 
measure or not. We can start with an order k > d/2. If (|4.1I) is infeasible, then 
we know y admits no -fT-measures, by Proposition 14.11 (ii). If (|4.1I) is feasible, we 
solve it for a minimizer w* if it exists. If w* has a flat truncation w*\2t with 
2t > deg(.A), then y and w*\2t commonly admit a finitely atomic i^-measure. If 
such a flat truncation does not exist, we increase the order k by one, and repeat 
the above. This produces the following algorithm for ,4-TKMPs. 

Algorithm 4.2. A semidefinite algorithm for ,4-TKMPs. 

Input: An „4-tms y, an even degree d > deg(A), a semialgebraic set K as in (11. ip 
and p > with K C B(0,p). 

Output: A finitely atomic if -representing measure for y, or an answer that it does 

not exist. 

Procedure: 

Step 0: Choose a generic R £ Ti n ^ and let k := d/2. 

Step 1: Solve (|4.ip . If (|4.ip is infeasible, output the answer that y admits 
no if-measures, and stop. If (|4.1|) is feasible, get a minimizer w* ,k . Let 
t := min{dK , deg(A)} . 
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Step 2: Let z := w*' \2t- Check whether the rank condition (|2.4|) is satisfied 

or not. If yes, go to Step 3; otherwise, go to Step 4. 
Step 3: Compute the finitely atomic measure y, admitted by z: 

[i = ci5{ui) + • • • + c r S(u r ), 

where r = rankM t (z), each m <E K and c$ > 0. Output /i, and stop. 
Step 4: lit < k, set t := t + 1 and go to Step 2; otherwise, set k := k + 1 and 
go to Step 1. 

For the input, we typically choose d = 2|~(deg(_4) + l)/2] , which is the minimum 
as required. In Step 0, the genericity means that R is chosen in S„ )t j\0, for a 
set C M N d having zero Lebsgue measure. In implementations, we choose R = 
[x]^ 2 G T G[x]d/2 with G a random square matrix obeying Gaussian distribution. 
In Step 2, the rank condition (|2.4[) is usually checked by using numerical ranks, 
due to computer round-off errors. In the numerical experiments of this paper, we 
evaluate the rank of a matrix as the number of its singular values that are greater 
than or equal to 10~ 6 . In Step 3, the method in [21] can be used to get a r-atomic 
if- measure for z. Algorithm 14.21 can be easily implemented by using the syntax in 
software GloptiPoly 3 [35] and SeDuMi [37]. Example 14.41 shows how to do this. 



Remark 4.3. If y admits a if-measure, Algorithm 14 . 2 1 can produce a r-atomic K- 
representing measure with r < \A\, for almost all R € either asymptotically 

or in finitely many steps (cf. Section [5j. The obtained r may not be minimum. It is 
typically quite difficult to find a representing measure whose support is minimum. 
This is an important future work, and we do not focus on it here. However, this 
problem can be treated in some way. For instance, Algorithm 14.21 can be applied 
repeatedly for a certain number of times, say, N. In each time, a different R is 
generated, and we typically get different r-atomic measures. Among these N times, 
we choose the measure whose support is minimum. Our numerical experiments 
show that this often produces a r-atomic measure with r equal or close to the 
minimum. Of course, this is heuristic, and there is no theoretical guaranty. 

Example 4.4. Consider the yl-tms y with (a,y a ) given as: 



Vo 



1/3 



1/3 











1/9 



1/15 



1/15 



This y admits a measure supported on the square [- 
average of x a on [—1, l] 2 . Clearly, h = $, g = (1— x\, 1— 



1,1] 2 , since each y a is the 
r?,)and [-1.1] 2 C B(0,V2). 
In (]4.ip . the tuple gs can be replaced by g, because ^k(gB) = ^k(g)- We apply 
Algorithm [12] to this ATKMP. The order k = 4 is typically enough to get a flat 
extension. This can be done by the syntax in GloptiPoly 3 [22] and SeDuMi [37] 
as follows: 

mpol x 2; 

Amon = [x(l)~2 x(2)~2 x(l)~2*x(2) x(l)*x(2)~2 ... 

x(l)~2*x(2)~2 x(l)~4*x(2)~2 x(l) ~2*x(2) "4] ; 

y=[ 1/3 1/3 1/9 1/15 1/15]; 

conmom = [mom (Amon) ==y] ; 

K = [l-x(l)~2>=0, l-x(2)~2>=0] ; 

bracx = mmon(x,0,4); G = randn (length (bracx) ) ; 
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R = bracx'*(G'*G)*bracx; k = 4; 
P = msdp(min(mom(R) ) ,K,conmom,k) ; 

[A,b,c,S] = msedumi (P) ; 

[xsol,ysol,info] = sedumi(A,b,c,S) ; 
dvar = c-A'*ysol; 

Mw = mat( dvar(S.f+l:S.f+S.s(l)~2) ); 

In the above, Mw is the moment matrix Mk{w). By Remark 14.31 wc run Algo- 
rithm 14.21 for a couple of times. In each time, the computed Mw satisfies the rank 
condition (|2.4[) . and we got a r-atomic measure with r < \A\ = 7, by the method 
in [21j . The smallest r we got is 3, which occurred in the representing measure 
Y^=i CiS(ui) with Ui and Cj given asQ 



Ui 


(-0.8524, 0.8910) 


(0.2109, 0.8873) 


(0.6697, -0.3743) 


Ci 


0.1231 


0.2061 


0.5233 



Example 4.5. Consider the .4-tms y with (a,y a ) given as: 



1/5 



1/15 



1/15 



3/105 



1/105 



1/7 



The above y admits a measure on § 2 , because each y a is the average of x a on S 2 . 
The sphere § 2 is defined by h = (||x||| - 1) and g = 0. Clearly, S 2 C B(0, 1). In 
(|4.ip . the tuple gs = (1— |M||) can be replaced by g — 0, because ^ k{gB)^Ek{h) = 
^k(g) H Ek(h). Like in Example 14.41 we run Algorithm 14.21 for a couple of times. 
In each time, we got a r-atomic measure with r < |^4| = 6. The smallest r we got 
is 2, occurring in the representing measure J2i=i c i$( u i) with Ui,Ci given as 





(0.3434 0.4542 0.8221) 


(0.8999 0.2550 -0.3539) 


Ci 


0.4610 


0.2952 



Example 4.6. (random instances) We apply Algorithm l4.2l to solve some randomly 
generated ATKMPs. Let K = 5(0,1) be the unit ball in R n . For each pair 
(n, m) from {(2, 10), (3, 8), (4, 6), (5, 4)}, we randomly generate a subset A C 
with cardinality 10, 20, 30 respectively. For each A, generate N := (™^ m ) points 
randomly from the ball B(0, 1), say, u\, . . . , un, and let y — ci[ui]^ + - • - + cn[un]a> 
with Ci > random. This is because if a tms in ]R N ™ admits a measure, then it 
must admit an A-atomic measure (cf. [1]). 



(n, m) 


\A\ = 10 


\A\ = 20 


|^| = 30 


(2, 10) 


4,5,6,7,8 


8,9,10,11,12,13 


9,10,11,12,13,14,15 


(3, 8) 


4,5,6,7,8,9 


9,10,11,12,13,14,15 


11,12,13,14,15,16,17,18 


(4, 6) 


5,6,7,8,9 


8,9,10,11,12,13,14 


10,11,12,13,14,15,16,17 


(5, 4) 


3,4,5,6,7,8 


7,8,9,10,11,12 


10,11,12,13,14,15,16 



For each triple (n, m, |-4|), we generate 100 random instances. For every generated 
instance, Algorithm 14.21 returned a r-atomic representing measure. The values of 
obtained r are listed in the above table. They are all smaller than \A\. This is 
justified by Proposition 15. 2f ii). 



Throughout the paper, only four decimal digits are shown for supports and weights. 
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5. Convergence Analysis 

In this section, we analyze the convergence of Algorithm 14.21 Two kinds of 
convergence will be investigated: asymptotic convergence, and finite convergence. 
For asymptotic convergence, we mean that there exists t such that the truncated 
sequence {u>*' fc |2t} (w*' k is a minimizer of (|4.1|) with order k) is bounded and all 
its accumulation points are flat extensions of y, if y admits a X-measure. For finite 
convergence, we mean that there exists k such that, either ()4.1[) is infeasible, or there 
exists t such that the truncation w* ,k \2t is flat. We begin with some properties of 
the semidefinite relaxations (|4.1[) - (|4.2I) . 

5.1. Properties of the semidefinite relaxations. 

Proposition 5.1. Let d > deg(_4) be even, and K C B(0,p) be as in U.l)) . Suppose 
y G M. A admits a K -measure. 

(i) If R lies in the interior ofYi n _d, then, for all k > d/2, is feasible and 
has a minimizer, and jp - JT^P have the same optimal value. 

(ii) Let Rmin be the optimal value of 13.3]) . Then, the optimal values of 
and are smaller than or equal to Rmin> for all k > d/2. 

(iii) Suppose R lies in the interior o/£„ i£ /. Then, there exists a constant C ~ 
C{R) such that, for all w that is a minimizer of lty with order k, 

(5.1) IMaHa < {1+ p 1 + ■■■+ P 2t )C, t = 0,l,...,fc. 

Proof, (i) Let \x be a finitely atomic i^-representing measure for y, which must 
exist by Proposition 13.31 Then the tms / R -[x] 2 fc d A i is feasible for (|4.1I) . So, (|4.ip is 
feasible for all k > d/2. Since R € mi(S nj d), for allp G ^[x]a with tiny coefficients, 
R — p G £„ d C l2k(h) + Qkigs)- This means that the zero polynomial is an 
interior point of (14. 2[) . Hence, strong duality holds, i.e., (|4.1[) has a minimizer, and 
(|4.ip - (14.2j) has the same optimum value, for all k > d/2 (cf. [2 Section 2.4]). 

(ii) Since y admits a A"-measure, the feasible set of (|3.3j) is nonempty. Let z be 
an arbitrary feasible point of (|3.3|) . Then z admits a finitely atomic K- measure /x, 
i.e., z = J K [x]ddiJ,. Let w = J K [x}2kd.ti- Clearly, w\a — y and is feasible for (|4.1|) . 
and (R,z) = (R,w). That is, (|4.1j) is a relaxation of (|3.3|) . So, the optimal value 
of (|4.1[) is at most R m in- This is also true for the optimal value of (|4.2[) , since it is 
not bigger than that of (|4.1I) . 

(iii) Since R G mt(E nj d), i? — e G for some e > 0. So, we have 

< (R - e,w) = (R, w) - e(l, w). 
(Cf. |31[ Lemma 2.5].) Since w is a minimizer of (|4.1j) . item (ii) implies that 

w = (l,w) < (R,w)/e < R m in/e- 
The membership w G $k(9B) implies L^\w) >z (b = p 2 — ||a;|||). So, 

p^dlxll^- 2 ) - & w {\\x\\f) > 0, i = 0, 1, . . . , k. 
(Cf. [311 Lemma 2.5].) A repeated application of the above gives 

^ w (\\x\\f) < p 2t w , i = 0,l,...,fc. 
Since Mk{w) > 0, for each t = 0, 1, . . . , k, 

t 

IM2J2 < ||Af t («;)|| F < Trace(M t (w)) = 

i=0 |a|=i 
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53 w 2a = ^ W (J2 x2a ) < ^(Mf) < p 2t wo- 

\a\=i 

Let C = Rmin/e, then the inequality (|5.ip holds. □ 

Proposition 5.2. Let K C H." be compact and d > deg(„4) be even. Suppose 
y £ R- 4 admits a K -measure. If R is generic in £ n ,d, then we have: 

(i) The problem i3.3\) has a unique minimizer. 

(ii) //, for some k, a minimizer w*' k of |^.-?| ) has a fiat truncation w* ,k \ 2t 
(2t > d), then the measure admitted by w* ,k \ 2 t * s r-atomic with r < \A\. 

Proof. The boundary of £ n ,d has zero Lebsgue measure in the space R N S. It is 
enough to prove the items (i) and (ii) if R is generic in the interior of £ n ,d- Let 

Se := {R £ S n , d : R-l/ee E M , \\R\\ 2 <£}, £=1,2,.... 

Clearly, int(E nid ) = [J e>1 int(Si). It is sufficient to prove that for each I = 
1,2,3,..., if R is generic in Se then the items (i) and (ii) are true. 

(i) For every R £ Se, we have l/£ < R(x) < £\\ [x] d \\ 2 for all x £ R™. Choose a 
finitely atomic measure ^* £ meas(y, K). Let Mi := J x £|| [x] d 1 1 2 d^*, a constant 
independent of i?. Clearly, J„[a;]dd^* is feasible in (I3.3p . and 

#mi„ < / i?di^* < Mi. 
J K 

For all i? € Si, (|3.3p has a minimizer z* (cf. Proposition 13. 6[) , and z* satisfies 

z* /l<(R,z*) = R min <M 1 . 

Note that z* is feasible for (|4.ip with order A; = rf/2. As in the proof of Proposi- 
tion [STTJiii), we can get 

||**||a <M 2 := (1 + p 2 + ■ ■ ■ + p d )lM x . 

Thus, for all R £ Se, (j3.3p is equivalent to 

(5.2) min (R,z) s.t. z\ A = y, \\z\\ 2 < M 2 , z £ @ d {K). 

z 

The feasible set of (|5.2p . denoted by F, is a nonempty compact convex set. A linear 
functional (R, z) has a unique minimizer on F if and only if R is a regular normal 
vector of F, or equivalently, (R, z) has more than one minimizer on F if and only 
if R is a singular normal vector of F (cf. Schneider (36J Section 2.2]). Let O be 
the set of all singular normal vectors of F. Then 9 has zero Lebsgue measure in 
the space R N S (cf. Schneider [36l Theorem 2.2.4]). If R £ S e \Q, then (JOJ), as well 
as p.3p . has a unique minimizer. So, if R is generic in Se, then (|3.3I) has a unique 
minimizer. 

(ii) Since w*' k \ 2t is flat, it admits a finitely atomic if-measure, and so does 
w*' k \ d when It > d. Thus, w*- fe | d is feasible in (|3~3|) . and (R, w*' k \ d ) > R mm . By 
Proposition I5.1f ii). we know (R, w*' k \ d ) = (R,w*> k ) < R rnin . Hence, (R, w*' k \d) = 
Rmin and w*\ d is a minimizer of (|3.3p . By item (i), we know, if R is generic in 
S nj d then (|3.3p has a unique minimizer. Therefore, for a generic R £ £ n ,d, w*\ d is 
the unique minimizer of f|3.3|) and is an extreme point of £ d (y, K). By Lemma l3.5[ 
every measure admitted by w* ,k \ d must be r-atomic with r < \A\. Clearly, every 
measure admitted by w* ,k \ 2 t is also admitted by w*' k \ d , and thus must be r-atomic 
with r< |^| . □ 
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5.2. Asymptotic convergence. Our main result in this subsection is: 

Theorem 5.3. Let d > deg(_4) be even and K C B(0,p) be as in U.l)) . Suppose 
y £ R" 4 admits a K-measure. In Algorithm ^, if R is generic in X n> <i, then J^.-?| ) 
has an optimizer w*' for every k > d/2 and we have: 

(i) For all t big enough, the sequence {w*' k \2t} is bounded and all its accumu- 
lation points are flat. Moreover, each of the accumulation points admits a 
r -atomic measure with r < \A\. 

(ii) In item (i), if, in addition, d is also big enough, then the sequence {w* ,k \2t} 
converges to a flat tms. 

Proof. Since R is generic in E ni< j, we can assume R £ mf(E„ l( j). By Proposi- 
tion l5.1f i). (|4.ip has an optimizer w*' k for every k > d/2. 

(i) By Proposition 15. If iii) ■ there is a constant C = C(R) such that 

\\w*' k \2th<(l + p 2 + --- + p 2t )C for all i = 0,l,...,fc. 

So, the sequence {iu*' fe |2t} is bounded, for any given t. Generally, we can assume 
p < 1 (otherwise, scale K to get p < 1). Let uj be an accumulation point of 
{w*' k \2t}- We can generally further assume w*' k \2t — > v as k — > oo. 

First, we prove that the truncation u)\d is a minimizcr of (pOjl . Each w*< k could 
be thought of as a vector in R N °° by adding zero entries to the tailing. The set R N °° 
is a Hilbert space, equipped with the inner product 

(u, v) :— u a v a , Vit, wSR N °°. 

a 

In the above, we know that for all k 

\\w*' k h = IK'IklU < (i + p 2 + ■ ■ ■ + P 2k )c < c/(i - p). 

The sequence {w*' k } is bounded in R N =°. By Alaoglu's Theorem (cf. TJ Theo- 
rem V.3.1] or [25J Theorem C.18]), it has a subsequence {w* ,kj } that is convergent 
in the weak-* topology. That is, there exists w* £ R N °° such that 

(c,w*' kj ) -t (c,w*) as j -t oo 

for all c £ R N =°. If we choose c as (c, w) = w a , for each a, then 

(5.3) w*^\ a ~>w*\ a . 

Since w*' k \ 2 t — > w, the above implies w*\ 2 t = w. Note that w*' k] £ {gB)^E kj (h) 
for all kj. For each r = 1, 2, . . ., if kj > 2r, then 

4?( w(fcj) ) = (1 < * < mi), L^(w^) hO(0<i< ma). 
Hence, (|5.3j) implies that for all r = 1, 2, . . . 

L h?(™*) = 0(1 <*< mi), ^ 0(0 < * < n»a)- 

This means that w* £ R N 'i is a full moment sequence whose localizing matrices 
of all orders are positive semidcfinite. By Lemma 3.2 of Putinar [25], w* admits a 
if-measure. So, the truncation oj\d = w*\d is feasible for (|3.3[) . and 

Rmin < (R,W*\d) = (#,™*) = 

By Proposition 15 . 1 f ii) . we know (R,w* ,k ) < R m in for all k. Thus, 

(R,w*) = lira (R,w*< k J) < R mm . 
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Hence, {R,u>\d) = Rmin, and uj\ d is a minimizer of f|3 . 3[) . 

Second, we prove that if t > \A\cIk then the truncation oj\2t is flat. By Propo- 
sition [521 if R is generic in £„ j( j, then (|3.3|) has a unique minimizer, which must 
be u\d, by the above. So, w|d is an extreme point of £ d {y,K) which is the feasible 
set of (|3.3p . Let p* be a If- representing measure for w, which must exist because 
u) = w*\2t- Then p* € meas(a>|<i,if). By Lemma 13. 5[ p* must be finitely atomic 
and |supp(p*)| < \A\. Note that 

rankA/ (w) < rankMd K (u;) < • • • < va,nkMd K \j.\( u} ) < ' ' ' > 

rankMi(w) < |supp(^*)| < |.4|, i = 1,2, .... 

There must exist i < \A\ such that 

rankM(^ 1)dK (w) = rankM Mf( (w). 

So, the truncation uj \ 2id K is flat. By Theorem l2.2[ it admits a unique if- representing 
measure. Note that p* is a if -representing measure for u}\2td K , and it is unique. 
Clearly, uj is an extension of io\ie.d K i and every measure admitted by uj is also a 
representing measure for uj\2id K - Therefore, p* is also the unique if- representing 
measure for u). Clearly, p* is r-atomic and r — rankil^dK^)- 

|suppO*)| =rankM^ K (cj) < 
rank Mt dK+1 (u) <■■■< rankAi t (w) < |supp(/x*)|, 

and we have rank M t -d K (w) = rankAi t (a;), i.e., uj is flat. 

Third, in the above, we have indeed shown that if fj,* is a if -representing measure 
for any accumulation point oj, then n* is r-atomic with r < \A\. 

(ii) It is enough to show that if t > d/2 > \A\dx then {w*' k \2t} has a unique 
accumulation point. We continue the proof of the item (i). By Proposition 15.21 
if R is generic in S„ ; d, then (|3 . 3[) has a unique minimizer, say, z*. Let uj be an 
arbitrary accumulation point of {w*' k \2t}- We have shown that cj\d is a minimizer 
of (|3.3|) in (i). So, cj| d = z*. If d > 2\A\d K , we have shown that z* is flat in (i). By 
Theorem l2.2[ z* admits a unique if -representing measure, say, v* . Since uj\d — z*, 
every if- measure admitted by uj must also be a if -representing measure for z*. 
Hence, v* must be the unique if -representing measure for uj, i.e., uj — f^x^tdv* ■ 
This shows that {w*' k \2t} has a unique accumulation point, which is f„ [x]2td;/* . □ 

Algorithm 14 . 21 is guaranteed to converge with probability one, for all ^4-tms y that 
admits a if-measure. If, accidently, a bad R is generated such that Algorithm 14.21 
fails to converge, we can choose a different generic R € S n;( j- Indeed, this never 
happened in our numerical experiments. Theorem 15.31 guarantees that we almost 
always succeed by choosing R from Ti n d- 

5.3. Finite convergence. Now we characterize when Algorithm l4.2l has finite con- 
vergence. Denote g„ l2+ i := p - \\x\\%, then g B = (gi, ■ ■ ■ , g m2 , 5m 2 +i)- For an index 
set J — {ji, . . . , j r }, denote gj := (g n ,...,g jr ). 

For a polynomial /, denote by f hom the homogeneous part of / with the highest 
degree. If p = (pi, . . . ,p r ) is tuple, denote p hom := (p\ om , . . . ,p^ om ). We denote 
by Jac(p)\ u the Jacobian of p evaluated at the point u. The discriminant of p hom : 
denoted as A(p hom ), is a polynomial in the coefficients oip hom , such that A(p hom ) = 
if and only if p hom (x) = has a nonzero solution «gC° with rank Jac(h hom )\ u < 
r. We refer to [30l Section 3] for discriminants. 
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Assumption 5.4. For any J C [m,2 + 1] with Vk(h,gj) ^ %, the coefficients of h 
and g B satisfy A(h hom , g h j om ) ^ 0. 

Assumption 15.41 holds generically, because it requires (ji hom ^g hom ^ to lie in an 
open dense set. The finite convergence of Algorithm 14.21 is as follows. 

Theorem 5.5. Let K C B(0,p) be as in U.l)) . and y be an A-tms. Let w*' k be an 

optimizer of j4-l\ l with order k, if it exists. 

(i) Suppose M\x\ A is K-full. Ify admits no K-measures, then f^. lty is infeasible 
for all k big enough. 

(ii) Suppose measly, K) ^ and Assumption \5.4\ holds. If \3.J$ has a maxi- 
mizer p* with R — p* G iafei (h) + Qk t (gB) for some k\, and if R is generic 
in Tj n: d 7 then w*' k \2t is flat for some t > d/2, for k big enough. 

(iii) Suppose achieves its maximum for k big enough. If the truncation 
w*' k2 \2t is flat for some k 2 and t > d/2, and if R G int(Yj n ^d), then \3.$ 
has a maximizer p* with R — p* G Iiks Q 1 ) + Qk 3 (gB) for some k^. 

Proof, (i) If R[x] A is K-full and y & M A (K), by Corollary EU there exists p 1 G 
R[x]^ such that Pi\k > 0, (pi,y) < 0. By the K- fullness of M. [x] A , there exists p 2 G 
R[a;]yl with p 2 \x > 0. So, for e > small, p :— p\ + ep 2 satisfies p\x > 0, (p, y) < 0. 
Let r]o > be such that (R + T]oP)\k > 0. By Theorem 12. 4[ both R + rj p and p 
belong to G htrih) + QtAgB), for some t\. Hence, for all rj > tjq and k > t%, —rjp 
is feasible for (I4.2[) . because 

R + rjp = R + r/op + (rj - r/ Q )p G hk(h) + Qk(g B )- 

Note that (—r]p,y) — > +00 as r\ — > +00. So, the optimal value of (|4.2[) is +00. By 
weak duality, its dual problem (|4.1j) must be infeasible for all k > t\. 

(ii) By Proposition 15. lf i). (|4. 1 [) and (|4.2[) have the same optimal value, for all k, 
if R G int(Ti n ^d) (this is true if R is generic in £ ni <j). For all k > kx, p* is feasible 
for ipnjl . So, if k > h, then (R,w*- k ) = (p*,y) = (p*,w*- k ), and (R-p*,w*< k ) = 0. 
Let q:=R-p*. Then 

(q,w*' k ) = Wk>ki. 
Clearly, q is nonnegative on K. Let z* be a minimizer of (|3.3[) and /1 G meas(z, K). 
Then, by Proposition 13.61 (i) (note R\k > if R is generic in £„.<j), 

0=(R,z*)-(ja*,y) = (q,z*) = f qdp. 

Thus, q vanishes on supp(/i), and has a zero on K. 
Consider the optimization problem: 

(5.4) min q(x) s.t. h(x) = 0, ffsfe) > 0. 

X 

The fe-th order SOS relaxation for (|5.4p is 

(5.5) 7 fc := max 7 s.t. q - 7 G hk(h) + QkigB)- 
Its dual problem is 

(5.6) min (q,w) s.t. w G <&/c(5b) H E k (h),w a = 1. 

The minimum of g over if is 0, and 7fe = for all k > k\. Thus, the sequence 
{7fc} has finite convergence. The SOS program (|5.5|) achieves its optimal value for 
k> ki, because q G 7 2 fe! (ft) + Qkt (gB)- 
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Since d > deg(^), (R - p*) horn = R hom . Under Assumption El for every 
J C [m 2 + 1] with Vk(h,gj) ^ 0, A(f,h hom ,g>} om ) is not constantly zero in / 
(cf. [IB Theorem 3.2]). So, if R is generic in E„ jd , then A(R hom , h hom , g'} om ) ^ 
for all such J. By Proposition IA.1I in the Appendix, (|5.4p has only finitely many 
critical varieties and Assumption 2.1 in [31] for (|5.4j) is satisfied^. If = 0, then 
w*' fc = for all a because M k (w*' k ) >z 0, and w*' k is clearly flat. If WQ k > 0, 
we can scale w*' k so that Wq' = 1. Then w* ,fc is a minimizer of (I5.6[) because 
(q,w*> k ) = for all fc > fej. By Theorem 2.2 of 31j, w*> k has a flat truncation 
w*' k \ 2 t if & is big enough. Indeed, w*' k \ 2 k-2 is hat (cf. Remark 2.3 of [H]). So, 
there is a flat truncation w* ,k \ 2t with £ > d/2 if fc is big enough. 

(iii) Suppose w*' k2 \ 2t is flat and 2t > d. Let R m in be the optimal value of (|3.3I) . 
then by Proposition 15 . If ii) . i? m m > (R,w* ,k2 ). On the other hand, the truncation 
w* M \ d is feasible in (|33| . so i? mi „ < (R,w* M ). Thus, i? mm = (R,w* M ). Indeed, 
we have i? m in = (R 7 w*' k ) for all k > k 2 . By assumption, (|4.2j) has a maximizer p* 
at a big order, say, > fea. Then R—p* e I 2 k 3 (h) + Qk 3 (gB)- Clearly, p* is feasible 
for (|3.4| because (Ji— > 0. By Proposition l5.1f i). the optimal value of (|4.2j) is 
also equal to R m in, if -R G int(S„ iC j). So, (p*,y) — (R, w* ,fc3 ) = R m in- The optimal 
values of (|3.3[) and (|3.4I) are equal, by Proposition [3^1) if i? € int(Y> n ^)- Hence, 
p* is a maximizer of Q3.4p . □ 

Remark 5.6. By item (i) of Theorem 15. 5[ if is if-full and y admits no K- 

measures, then (|4.1I) is infeasible for some k, for any R (we don't need R £ E nc j). 
When is not if -full and y admits no if -measures, it is not clear whether 

or not there exists k such that (|4.ip is infeasible. This is because there does not 
exist a characterization like Theorem 13. II for the membership in 3% a (K) if K[x]^ is 
not if -full, to the best of the author's knowledge. Theorem 13.11 and Corollary 13.21 
assume if-fullness of R[x]^t. In many applications, M[x]_4 is often if -full. On the 
other hand, if y admits a if-measure, no matter R^]^ is A'-full or not, for a generic 
R G Sn,d, Algorithm 14.21 will find a finitely atomic if -representing measure for y, 
cither asymptotically or in finitely many steps. 

Remark 5.7. a) When int(K) ^ 0, (|4.1j) has interior points, and thus (14. 2p 
achieves its optimal value, for every order k (cf. [20j [25|). b) By Theorem 15.51 
(ii) and (iii), the condition R — p* 6 1(h) + Q(gs) is almost necessary and sufficient 
for finite convergence to occur, modulo some general technical assumptions, c) If 
a polynomial / is nonnegative on if, then / G 1(h) + Q(.gs), under some general 
conditions (cf. [32]). So, the condition R — p* G 1(h) + Q(gs) is often satisfied. 
Thus, it is very likely that Algorithm 14.21 has finite convergence. Indeed, the finite 
convergence occurred in all our numerical experiments. 

6. Applications 

In this section, we show how Algorithm 14.21 can be applied to solve CP/SOEP- 
decomposition problems and the standard truncated if- moment problems. 



2 In |31| , polynomial optimization problems with only inequality constraints were discussed. 
If there are equality constraints, Assumption 2.1 in I31| can be naturally modified to include all 
constraining equations, and Theorem 2.2 of | 31| is still true, with the same proof. 
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6.1. Completely positive matrices. Recall that a matrix C £ S n is completely 
positive if there exist iti, . . . ,u r £ R" such that 

(6.1) C = uiu( + • • • + U r Ur ■ 

If (|6.ip holds, we say C is a CP-matrix. The number r is called the length of the 
decomposition (I6.1[) . The smallest r, for which (|6.1[) holds, is called the CP-rank of 
C (cf. [S]). Let Cp(n) be the cone of n x n CP-matrices. Clearly, C £ Cp(n) if and 
only if C = BB T for a nonnegative matrix B (i.e., every entry of B is nonnegative) . 
So, every CP-matrix must be positive semidefinite, but typically not vice versa. 
The dual cone of Cp(n) is the set ofnxn copositive matrices (a matrix A £ S n is 

copositive if x T Ax > for all x £ R"). Let A„ = {x £ R™ : xi H h x„ = 1} be 

the standard simplex in R™. 

Completely positive and copositive matrices have wide applications in optimiza- 
tion, like approximating stability numbers (cf. |13j ) or solving nonconvex quadratic 
programs (cf. [B]). Checking the membership in Cp(n) is NP-hard (cf. [14 ]). We 
refer to the survey |16] by Diir and the book [5] by Berman and Shaked-Monderer. 
If C is not a CP-matrix, how can we get a certificate for that? If it is, how can we 
get a CP-decomposition for CI When C is acyclic or circular, Dickinson and Diir 
|15) showed that checking complete positivity can be done in linear-time. In the 
general case, Berman and Rothblum [4 showed that checking complete positivity 
and computing CP-ranks can be done by using Renegar's algorithm on quantifier 
elimination |33j . This method is symbolic, typically runs in exponential time, and 
is usually very expensive to implement. In the prior existing work, there are no 
much efficient numerical methods for solving general CP-decomposition problems, 
in the author's best knowledge. 

Clearly, C is a CP-matrix if and only if 

(6.2) C = QiUiuJ + • • • + Q r ll r Ur , 

for some u\,...,u r £ A„, Qi,...,Q r > 0. Every symmetric matrix C can be 
identified by the vector consisting of its entries 

Let Q n = {a £ N n : \a\ = 2}. Then c is a Q„-tms, and (16.21) is equivalent to 

C = Ql{ui\Q n H h Qr[u r ]Q n - 

Clearly, if C £ Cp(n), then c admits a A n -measure. Conversely, if c admits 
a A„-measure, then c also admits a finitely atomic A n -measure like the above 
(cf. Proposition I3.3j) . and C £ Cp(n). Thus, the CP-decomposition problem is es- 
sentially an ,4-TKMP with A = Q n , K = A n . The simplex A„ is in the form ([L~T]) . 
with h — (l T x — 1) and g = (x%, . . . ,x n ) (1 denotes the vector of all ones), and 
A„ C B(0, 1). Note that R[x] c „ is A„-full. 

By the above, the CP-decomposition problem can be solved by Algorithm 14.21 
If C $ Cp(n), then Algorithm 14.21 will return a certificate for this (i.e., (|4.1[) is 
infeasible for some k), by Theorem [53^1). If C £ Cp(n), then we can asymptotically 
get a flat extension of c, for almost all R £ S n> d (d > 2 is even), by Theorem 15.31 
Moreover, we can likely get it within finitely many steps (cf. Remark |5 . T|) . Indeed, 
finite convergence occurred in all our numerical experiments. After getting a flat 
extension of c, we can get a r-atomic A„-representing measure for c, which then 
produces a CP-decomposition for C. 
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Example 6.1. Consider the matrix: 



C 



6 4 12 
4 6 1 
10 3 1 
2 112 
2 3 2 1 



We apply Algorithm 14.21 to the corresponding Qs-tms 

c:= (6, 4, 1,2, 2, 6, 0,1, 3, 3, 1,2,2, 1,5). 



To get a decomposition of small length, we run Algorithm 14 . 2 1 for a couple of times 
(cf. Remark I4.3|) . In each time, we got a CP-decomposition for C. The smallest 
length we got is 5, which occurred in the factorization C = BB T with 



B 



1.0911 2.0836 

0.0000 1.6456 

0.1488 0.0000 

1.0797 0.3606 

0.5830 0.0000 



0.0000 
0.0000 
1.0379 
0.8087 
0.0000 



0.3148 0.6076' 

1.8143 0.0000 

0.0000 1.3786 

0.2241 0.0000 

1.6535 1.3878 



The CP-rank of C is 5, because 5 = rankC < CP-rank C < 5. 
Example 6.2. Consider the matrix (cf. [5, Example 2.9]): 



C 



110 1 
12 10 
12 10 
12 1 
10 16 



It is positive semidefinite, but not completely positive (cf. [5]). We apply Algo- 
rithm 03] to verify this fact. It terminates at Step 2 with k = 2, because (|4.ip is 
infeasible. This confirms that C is not a CP-matrix. 



Example 6.3. (random instances) We apply Algorithm d21 to randomly generated 
CP matrices. If C G Cp(n), then C admits a CP-decomposition (|6.2j) with length 



r < \n(n + 1), by Caratheodory's Theorem. Indeed, it can be slightly sharpened 
to r < ±rankC*(rankC* + 1) - 1, if rankC > 2 (cf. [55]). Clearly, we always 
have r > rank(C). So, if C G Cp(n) and C has full rank, then n < r < cp(n) := 
\n(n + 1) — 1, for n > 1. For n = 2, 3, . . . , 8, we generate 50 instances, except for 
n = 8 (only 20 instances are generated). For each instance, generate N := |n(n+l) 
points randomly from A„, say, u%, . . . , ujy, and let C = ciUiu[ + ■ ■ ■ + cnUnu^ 



with Ci > random. For each C, we apply Algorithm 14 . 2 1 ten times and let r be the 
smallest length that is obtained. Algorithm 14.21 is able to get a CP-decomposition 
for all generated C. The obtained values of r are listed in the table: 



n 


2 


3 


4 


5 


6 


7 


8 


cp(n) 


2 


5 


9 


14 


20 


27 


35 


r 


2 


3 


4 


5,6 


6,7,8 


8,9,10 


11,12,13,14,15 



They are equal or close to the lower bound n (because rankC 
C), and is much less than the upper bound cp(n) for n > 4. 



n for generated 
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6.2. Sum of even powers (SOEP) of real linear forms. Recall that a form / 
of an even degree m is SOEP if for some real linear forms L\, . . . , L r 

(6.3) / = L? + ■ ■ ■ + L™ 

Let Q n ,m denote the set of all SOEP forms in n variables and of degree m. Reznick 
proved that Q n ,m is a convex cone with nonempty interior and its dual cone is 
the set of nonnegative forms in n variables and of degree m. We refer to Reznick 
[35] for SOEP forms. The number of sums, r, is called the length of the SOEP- 
decomposition (|6.3p . The minimum r for which (|6.3I) holds is called the width of 
/, and is denoted as w(f) (cf. |35j). The decomposition (|6.3j) is called minimum if 
r = w(f ). SOEP decompositions naturally have wide and interesting applications, 
like in Waring's problems, quadrature problems, sphere designs [35 . It is typically 
quite difficult to check whether a form is SOEP or not. As shown by Reznick [55] . 
when m > 4, a rational form / 6 Q n ,m may not have a decomposition (|6.3[) with 
all Li rational. Therefore, numerical methods are preferable in applications. In 
the prior existing work, there are no much efficient numerical methods for SOEP- 
decompositions, in the author's best knowledge. 

Denote := {a £ N n : \a\ = m}. We can write a form / of degree m as 



Denote / := (/ a ). So, / can be identified by the H™ -tms /. If / is SOEP and (IO 
holds, then we can write each Li = ^/Ci(ufx) with Ci > and 

Ui e §^ _1 := {x e S' 1 " 1 : l T x > 0}. 

Thus, we get 

r / 
f ^ ^ 



i=l q£H" 



The above is equivalent to the decomposition: 

(6.4) / = Ci [u x ]h» H hc r [u r ]H» ■ 

Clearly, if / is SOEP, then / admits a §"~ ^measure. Conversely, if / admits an 
§" _1 -measure, then / also admits a finitely atomic §" -1 -measure (cf. Proposi- 
tion [33]), and so / is SOEP. Hence, checking / £ Q n ,m is equivalent to determining 
whether the H^-tms / admits a S™ _1 -measure or not. The latter is a ,4-TKMP 
with A = H™ and K = S^ _1 . Note that R[x]u%, is §" _1 -full. 

SOEP-decomposition problems can be solved by Algorithm l4.2l The set S™" 1 is 
as in (HI]) with/i= (||a;||l-l) and g = (l T x). Clearly, S^ 1 C B(0,1). In g3]), the 
tuple gs can be replaced by g, because ^k(g)^Ek(h) = ^k{gB)<^Ek(h). If / is not 
SOEP, then / admits no S™ _1 -measures, and Algorithm 14.21 can give a certificate 
for this fi.e.. (|4.1l ) is infeasible for some order fc), by Theorem 15. 5f iV If / is SOEP, 
then we can asymptotically get a flat extension of /, for almost all R £ T, n ^ (d > m 
is even), by Theorem l5.3l Moreover, we can likely get it within finitely many steps 
(cf. Remark 15. 7p . Indeed, this occurred in all our numerical experiments. Once 
a flat extension of / is obtained, we can easily get an SOEP-decomposition for / 
from a finitely atomic representing measure for /. 
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X3 



+ 



X2+X3) + {xi + x 2 - x 3 T) ■ 
1/3, we can get an SOEP-decomposition of length 11 for qi/3, by the 



Example 6.4. Consider the sextic form 

qx := (xj + x\ + x\f - \{x\ + x\ + x%). 

It is SOEP if and only if A < 2/3 (cf. 35, p. 146]). For A = 2/3, by running 
Algorithm 14.21 a few times, we got an SOEP decomposition of length 10 for ^2/3: 

JS ((Xl + X 2 f + (X 2 + Z3) 6 + (Xl + X-sf + (Xl - X 2 f + (Xl - X Z f + (X 2 
gg ((ll + X 2 + X S ) 6 + (~X! + X 2 + X3) 6 + (ll 

For A 

same way. The lengths 10 and 11 are the smallest ones that we can get for g 2 /3 
and (71/3 respectively (cf. Remark 1 5. 7p . When A = 1, (|4.ip is infeasible for k — 4, 
and Algorithm 14.21 terminates at Step 2. This confirms qi $ Q?,fi- 

Example 6.5. (random instances) We apply Algorithm [472] to randomly generated 
SOEP forms. Let m be an even degree. If / € Q n ,m, then its width w(f) < 
N := by Caratheodory's Theorem. If / G int{Q n ^ m ), then its width 

w{f) > N := (™ + ^2 _1 ) ( cf - S3 Theorem 3.14(iv)]). If n = 2 or (n,m) = (3,4), 
then w(f) < Nq (cf. [35l Theorem 4.6]). So, for the above ranges of (n,m), we 
know the generic width is Nq. For other ranges of (n, m), if / is generic inside 
Qn,m : then Nq < w(f) < N. We consider (n,m) from the table: 



(n, m) 


(2,4) 


(2,6) 


(2,8) 


(2,10) 


(3,4) 


(4,4) 


(3,6) 


(3,8) 


gwidth 


3 


4 


5 


6 


6 


[10,35] 


[10,28] 


[15,45] 


r 


3 


4 


5 


6 


6,7 


12,13,14 


11,12 


17,18,19 



In the above, gwidth is No if n = 2 or (n, m) = (3, 4), and is the range [JVo, N] for 
other cases. For each pair (n, m) from the above table, we generate 50 instances. 
In each instance, generate points Ui, . . . ,un randomly from §™ _1 , and let / = 
Ci(uf x) m + ■ ■ ■ + C7v(u^x) m with Cj > random. For each generated /, we run 
Algorithm 14 . 2 1 for ten times, and choose r to be the smallest length of the obtained 
SOEP-decompositions. For all generated /, we got an SOEP-decomposition, and 
the values of obtained lengths r are listed in the above table. We can see that r is 
equal or close to the minimum. 

6.3. Standard truncated if-moment problems. When A = N 7 ^, the ^l-TKMP 
is specialized to the standard truncated if- moment problem (TKMP), which was 
originally studied by Curto and Fialkow [SJ |U [TTJ [T^] . Algorithm 14.21 can be nat- 
urally applied to solve TKMPs. The set R[a;] m is if- full, for any set K. We are 
interested in the case that K is compact. If a tins y e M. h "" admits no if -measures, 
Algorithm l4.2l can return a certificate for the nonexistence of representing measures; 
if y admits a if -measure, we can asymptotically get a flat extension of y, for almost 
all R € £„ d (d > m is even), by Theorems 15.31 moreover, we can likely get it 
within finitely many steps (cf. Theorem 15.51 and Remark |5.7[) . In the author's best 
knowledge, Algorithm 14.21 is the first numerical algorithm that can solve general 
TKMPs with a compact semialgebraic set if. 

Example 6.6. Consider the following tms in R N s; 

(1, 0, 0,1/3, 0,1/3, 0, 0, 0, 0, 1/5, 0, 1/9, 0, 1/5, 
0, 0, 0, 0, 0, 0,1/7, 0, 1/15, 0, 1/15, 0, 1/7). 
Its moments are listed in the graded lexicographical ordering. This tms admits a 
measure supported on [— 1, l] 2 , because its a-th moment is the mean value of x a 



on [— 1, 1] . We apply Algorithm 14.21 to this tms. In each time of running, we got 
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a r-atomic representing measure supported in [— 1,1] . After a repeated running, 
the smallest r we got is 10 (cf. Remark I4.3j) . which occurred in the representing 
measure J2i=i c i&i u i) with Ui and q given as: 





Hi Ci 


(-0.7983, -0.9666) 0.0318 


(-0.8710, -0.2228) 0.0947 


(-0.2155, -0.6211) 0.1698 


(-0.9175, 0.8643) 0.0328 


(0.5833, -0.8876) 0.0729 


(-0.4834, 0.4714) 0.1662 


(0.3541, 0.0676) 0.2054 


(0.9269, -0.3819) 0.0672 


(0.0841, 0.9294) 0.0717 


(0.8153, 0.6919) 0.0874 
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Appendix A. Generic Finiteness of critical varieties 
Consider the polynomial optimization problem 
min fix) 



(A 1 

{ ' ' \ s.t. hi(x) = (i € [mi]), gj(x) > (J e [m 2 + 1]). 

For each J C {1, . . . , rn 2 + 1}, denote 

Vj := {x G C™ : h(x) = 0,gj(x) = 0, isakJac(f,h,gj)\ x < m 1 + \J\}. 

The set of critical points of (|A. 1|) with the active set J is contained in Vj, which is 
called a critical variety of dHJ. We show that if the coefficients of f hom , h hom , g) om 
satisfy some discriminantal inequalities, then Vj is finite. We refer to t 30| Section 3]) 
for the definition of discriminants A. 

Proposition A.l. Let /, hi (i e [mi]),gj (J E [m,2 + 1]) € M.[x], and Vj be defined 
as above. For any J C [rri2 + 1], if 

Aif hom , h hom ,g h j om ) ^ 0, A(h hom ,g$ om ) ^ 0, 

then Vj is finite. 

Proof. Denote x :— (xo, x\, . . . , x n ) and by p the homogenization of a singleton or 
tuple of polynomials p. Let Uj be the projective variety in P™ (cf. [T^]) defined as 

(A.2) rank Jac(f,h,gJ)\ x < m x + \J\, h(x) = gj(x) = 0. 

Clearly, if u e Vj, then (l,u) G Uj. Suppose otherwise Vj is infinite, then Uj 
is positively dimensional. By Bezout's Theorem (cf. [19 ), Uj must intersect the 
hyperplane xq = in P™, i.e., (|A.2I) has a solution like (0, v) with =/= v £ C™. So, 
v is a solution to the homogeneous polynomial system 

(A.3) rank Jacif hom ,h hom , g h j om )\ x <m x + |J|, h hom ix) = g h j om ix) = 0. 

Since A(h hom , g h f m ) ^ 0, rank Jacih hom , g h j om )\ x = m x + \ J\ for all ^ x G 
Vc(/i /lom ,gj om ). The rank condition in CO]) implies that f hom iv) = (cf. [301 
Section 3]). Hence, w is a nonzero singular solution to 

f hom ix) = h hom ix) = g h j om ix) = 0, 

which contradicts A(/ hom , h hom ,g h j om ) ^ (cf. [30]). So, Vj must be finite. □ 
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